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Aubry-Andre-Harper (AAH) lattice models, characterized by reflection-asymmetric, sinusoidally 
varying nearest-neighbor tunneling profile, are well-known for their topological properties. We 
consider the fate of such models in the presence of balanced gain and loss potentials ±27 located at 
reflection-symmetric sites. We predict that these models have a finite VT breaking threshold only 
for speeifie loeations of the gain-loss potential, and uncover a hidden symmetry that is instrumental 
to the finite threshold strength. We also show that the topological edge-states remain robust in the 
T^T-symmetry broken phase. Our predictions substantially broaden the possible realizations of a 
PT-symmetric system. 


I. INTRODUCTION 

When is the spectrum of a non-Hermitian Hamilto¬ 
nian purely real or has only complex-conjugate pairs? 
Numerous authors have addressed this question, start¬ 
ing with Bender and coworkers who showed that (contin¬ 
uum) Hamiltonians invariant under combined parity- and 
time-reversal operations (T^T-symmetric) fit the bill [T]- 
[3]. Such Hamiltonians faithfully model open systems 
with balanced gain and loss, in which the parity opera¬ 
tor {V) exchanges the gain-region with the loss-region, 
whereas the time-reversal operator (T) transforms a 
gain-region into a lossy region. Concurrent with their ex¬ 
perimental realizations in coupled waveguides ms] , res¬ 
onators [ 9 ], microcavities m, and lasers EHni, discrete 
VT systems with a parity-symmetric tunneling term 
Hq = VHqV = Hq and a balanced gain-loss potential 
V = VTVVT 7 ^ have been intensely studied in the 
past five years p!5H2Q] . In particular, site-dependent tun¬ 
neling Hamiltonians, of interest for perfect-state-transfer 
and quantum computing [ 21 ] , have been theoretically [ 22 ] 
and experimentally [231125] explored. All of these experi¬ 
mentally investigated systems have been subject to a the 
stringent constraint of a reflection-symmetric tunneling 
amplitude profile. 

Generically, the spectrum of a T^T-symmetric Hamil¬ 
tonian H = Hq + U is real when the strength 7 of the 
balanced gain-loss potential is smaller than a positive 
threshold jpt set by Hq. The emergence of complex- 
conjugate eigenvalues at the exceptional point 7 = 7 pt 
is called T^T-symmetry breaking |26l|27|- When 7 > 
the eigenfunctions with complex eigenvalues become in¬ 
creasingly asymmetrical [28]. It has long been known 
that a purely real or complex-conjugate-pairs spectrum 
is equivalent to the existence of an antiunitary operator 
A = UT that commutes with the Hamiltonian H Ei- 
IM]. Thus, in principle, reflection-symmetry is not a 
necessary constraint, U ^ V. Indeed there are several 
proposals, based on supersymmetric quantum mechan¬ 
ics, for continuum models where the complex potential 
V{x) V ^ (~^) is not reflection symmetric (or antisym¬ 
metric) [32] [33]. Nonetheless, all experimental realiza¬ 


tions of T^T-symmetric systems to-date have abided by 
the reflection-symmetry constraint. 

On a separate front, one-dimensional Aubry-Andre 
Harper (AAH) models have been extensively explored in 
recent years in the context of their experimental realiza¬ 
tions in coupled waveguide arrays. Consider an Wsite 
tight-binding lattice with site-dependent tunneling pro¬ 
file tk = J [1 + Acos(27r/?/c + 0)]. Here J > 0 denotes 
the energy scale associated with the tunneling rate and 
dimensionless A characterizes the strength of the tunnel¬ 
ing modulation. When P = 1/2, this model is known 
as the dimer model or the Su-Schrieffer-Heeger (SSH) 
model; it describes the transport of charge carriers in 
acetylene [34| |35]. For a rational value of /?, the one¬ 
dimensional AAH model is related to the Hofstadter- 
butterfly problem [36] which describes the behavior of 
two-dimensional electron gas in a magnetic field in the 
presence of a periodic potential. When (3 is irrational, the 
AAH model describes one-dimensional quasicrystals [S3- 
ED]. For an infinite lattice, when /3 is rational, the tun¬ 
neling amplitude is periodic and the corresponding AAH 
model has robust topological edge states m and is re¬ 
lated to topological insulators [42ll45] . We emphasize 
that a lattiee with tunneling profile tk is not, in general, 
refieetion-symmetrie [46] . 

In this paper, we investigate the fate of Wsite AAH 
models - ones that are experimentally realizable in cou¬ 
pled waveguides or resonators - in the presence of one 
active, gain potential Pij at site tuq and a balanced loss 
potential —27 at its reflection-symmetric counterpart site 
mo = N 1 — mo @ 3 . Our four primary results are 
as follows, i) When P = 1/p, the model has a positive 
threshold 7 pt(^o) if and only if the lattice size N and 
the gain location mo both satisfy N 1 = 0 mod p and 
mo = 0 mod p. ii) When (3 = q/p is rational, where p, q 
are co-prime and g > 1 , the same pattern holds irrespec¬ 
tive of the value of g; when /3 is irrational, the threshold is 
zero, hi) When P = q/p, interspersed among its p bands, 
the model has p — 1 localized edge modes that continue 
to have real energies past the T^T-transition. iv) Our 
predictions are unaffected when the tight-binding lattice 
approximation is relaxed, and thus are valid in realis- 
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tic coupled optical waveguides. This work provides a 
pathway to investigate VT -symmetry breaking in lattice 
models with topological states. 

The paper is organized as follows. In Sec. [TT] we in¬ 
troduce the notation and summarize the properties of 
a finite, Hermitian AAH model. We present numer¬ 
ical results for the T^T-symmetry breaking threshold 
7 pt(^O: 0 ) for ^ wide range of lattice parameters, and 
summarize the findings. In Sec. |III[ we present a per¬ 
turbative analysis of the T^T-symmetry breaking thresh¬ 
old and show that, due to a hidden symmetry of the 
eigenfunctions of the AAH model, the threshold 7 pt is 
positive even when though the underlying system is re¬ 
flection asymmetric. In Sec. |IV| we consider the smallest 
such lattice, a dimer lattice with N = b sites. After 
an analytical solution, we present the dynamics obtained 
via beam-propagation method (BPM), which show that 
our predictions will remain valid in realistic samples. We 
conclude the paper with a brief discussion in Sec. |V[ 

II. LATTICE MODEL AND THE PT-PHASE 
DIAGRAM 

The Hermitian tunneling Hamiltonian for an A-site 
lattice with nearest neighbor tunneling and open bound¬ 
ary conditions is given by 

N-l 

Ho{X,l3, ^) = -Yl + 1| + + l)(fc|), (1) 

k=l 

t/e = J [1 + A cos(27r;d + (/))] , (2) 

where \k) denotes a single-particle state localized at site 
k. The parity (reflection) operator V on the lattice, in 
the site-basis, is given by Vab = ^ab where b = N -\-l — b. 
The time-reversal operator T = * where * denotes com¬ 
plex conjugation. The Hamiltonian is not, in 

general, invariant under the VT operation. The trivial 
exceptions are a uniform lattice, A = 0 , or a dimer model, 
(3 = 1 / 2 , with an even number of lattice sites. 

Since the tunneling function t/c(A, /3, (j)) is periodic in (3 
and 0 , without loss of generality, we consider P G [ 0 , 1 ) 
and (p G [0, 27r). It is also straightforward to show that 

Ho{-X,TP)=Ho{\TpT7r), (3) 

Ho{XA-TP)=Ho{X^T27r-cP). (4) 

Therefore, it is sufficient to restrict ourselves to A > 0 and 
P < 1 / 2 . The general band structure of the AAH Hamil¬ 
tonian i7o(A,yd,0) is highly intricate, where the number 
of bands is determined by T arid the locations of band 
degeneracies are determined by A and p [36H4T]. Note 
that when A > 1 , the tunneling amplitude tk changes sign 
from positive to negative along the lattice. In addition, 
for A > 1 and a rational p = g/p, the tunneling ampli¬ 
tude vanishes at /c = 0 mod p and p = arccos(— 1 /A). 
For such parameters, the A-site chain splits into pieces 
of size p, and the corresponding Hamiltonian Hq becomes 
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FIG. 1. (Color online): (a) Schematic of an AAH lattice with 
N — S sites, denoted by solid circles. The tunneling ampli¬ 
tude tk/J is periodic with period p = 1/0 = 3. Also shown 
are balanced gain-loss potentials ± 27 , denoted by red and 
blue solid circles respectively, at reflection-symmetric sites; 
the vertical dashed line is the lattice center, (b) Spectrum 
En{<p) of an A = 20 lattice with the same periodicity shows 
1/0 = 3 bands, each with [A0] =6 extended states; [x] de¬ 
notes the largest integer smaller than x. The two remaining 
mid-gap states are localized for all <p except 0 = { 7 r/ 3 , 47 r/ 3 } 
shown by dashed vertical lines. 


block-diagonal. In order to avoid such cases, whose be¬ 
havior can be trivially understood, we confine ourselves 
to modulation strengths 0 < A < 1. In the presence of 
reflection-symmetric gain-loss potentials ±27 at sites mo 
and mo, the lattice Hamiltonian becomes H = Hq + V 
with 

V = i7(|mo)(mo| - |mo)(mo|) = PTVVT^Vl (5) 

Figure [l] encapsulates the typical properties of Hamil- 
tonian Hq. Panel (a) shows the reflection-asymmetrical 
tunneling profile tk/J for an A = 8 site lattice with tun¬ 
neling modulation strength A = 0.5, inverse tunneling 
period 0 = 1/3, and phase 0 = 0. The neutral sites on 
the lattice are indicated by solid black circles, the solid 
red circle at mo = 3 denotes the gain site, and the loss site 
at its reflection-symmetric location mo = 6 is denoted by 
the solid blue circle. Panel (b) shows the energy spectrum 
En{p)/J for an A = 20 site AAH model with 0 = 1/3. 
In addition to the three bands of extended states that 
are expected at 0 = 1/3, there are two edge-localized 
states with energies that lie in the two band gaps. The 
midgap states are localized for all values of phases except 
0 = {7r/3,47r/3}, denoted by dotted vertical lines. 

These are generic features of the spectrum for 0 = 
g/p, which corresponds to the tunneling profile period of 
p, and lattice size A = Mp — 1. Each of the p bands 
has [A0] = (M — 1) extended states, and the remaining 
(p—1) midgap states are localized for almost all 0. When 
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+ 1 = 0 mod p, it is straightforward to show that 

VTH^{\ cl))VT = Ho{X, p, 27rP - 0 ). ( 6 ) 

Thus, Hq becomes VT symmetric if and only if 0 = 
{np^TTp + tt}. These are precisely the ^-values at which 
the midgap states become extended. In the following 
sections, we will see that the topological midgap states 
do not participate in the T^T-symmetry breaking and re¬ 
tain their localized character past the symmetry breaking 
transition. 




Gain-site 


FIG. 2. (Color online): Threshold jptIJ as a function of 
gain-potential location 1 < mo < N/2 and phase p for an A/" = 
59 lattice, (a) When P = 1/3, 7pt = 0 for all gain locations 
except mo = {3,6, • • • }. (b) When = 1/4, a positive 7 pt is 
obtained only when mo = {4, 8, • • • }. These results show that, 
contrary to naive expectations, a reflection-asymmetric AAH 
Hamiltonian has a positive threshold ^pT(jno^p)/J ~ 1 > 0. 

We now present the VT phase diagram for this model. 
Naively, the reflection-asymmetric nature of the tunnel¬ 
ing Hamiltonian Hq will imply, via perturbation theory, 
that an infinitesimal gain-loss potential, Eq.([^, will lead 
to a complex spectrum. This expectation is indeed con¬ 
firmed by numerical results for all lattice configurations 
except when /3 = g/p is rational, and the lattice size N 
and gain-site location mo satisfy A^ + 1 = 0 mod p and 
mo = 0 mod p. 

Figure]^ shows the VT threshold strength 7 pt (^0 5 P) 
for a lattice with A^ = 59 sites and tunneling modula¬ 
tion strength A = 0.5. Panel (a) shows the results for a 
tunneling profile with spatial period p = 3. The thresh¬ 
old strength is zero exeept when the gain loeation is an 
integer multiple of the tunneling period^ mo = 0 mod 3. 
Panel (b) shows that a similar behavior is obtained for 
P = 1/4. Note that this nonzero threshold results only 
for periods p such that N -\-1 = 0 mod p. Thus, for ex¬ 
ample, when p = q/p = q/7, the T^T-breaking threshold 
for an A^ = 59 site lattice is identically zero for any mo 
and any g > 1. In general, the nonzero threshold Jpt/J 
first decreases as the gain-potential site mo moves in from 


the end of the lattice and increases again as it approaches 
the lattice center, mo ^ N/2 [48]. These results are qual¬ 
itatively similar for large N and the maximum threshold 
strength remains the same in the thermodynamic limit, 
A^> 1. 

Figure shows the typical dependence of positive 
jPT on the tunneling period p and tunneling modula¬ 
tion strength A; in each case, only gain-potential loca¬ 
tions mo < N/2 that give rise to a positive VT thresh¬ 
old are considered. Panels (a) and (b) show the VT 
threshold in the (mo,0) plane for the same modulation 
strength A = 0.5 and lattice size N = 111. Consistent 
with the results in Fig. [^ the VT threshold varies non- 
monotonically with phase 0, and is generally maximum 
when the gain and loss locations are farthest apart or 
nearest to each other. As the tunneling period is in¬ 
creased from 1/P = 4, panel (a), to 1/P = 7, panel (b), 
we see that the region with appreciable threshold value 
shrinks in size but the maximum value of 7 pt does not 
alter substantially. 

Panel (c) in Fig. 1^ shows the dependence of the VT 
threshold 7 pt(^0 :^ on tunneling modulation A for a 
dimer lattice with N = 61. When p = 1/2, the tunnel¬ 
ing amplitude on adjacent bonds alternates between two 
values J(1 ^ A cos p) ISH |35], and the tunneling profile is 
not reflection-symmetric for an odd N. At A = 0.1, due 
to the small tunneling modulation, the threshold ypp is 
essentially independent of the phase 0, and its depen¬ 
dence on mo is similar to that for a uniform tunneling 
lattice; in particular, we see that ^pp/J 0.5 when 
the gain and loss sites are closest to each other [48]. As 
A increases, the VT threshold, which is proportional to 
the effective tunneling amplitude, is strongly suppressed 
when cos 0 = ±1, but remains unchanged from its A <C 1 
limit when \p\ ^ 7r/2 mod 2p. As an aside, we note for 
larger tunneling periods p > 2, the ^-dependence of the 
threshold ypp is not as easily characterized. Results in 
Figs. and might suggest that the threshold reaches a 
maximum nt p = {7r,7r + Tr^d}; however, that is not true 
for all modulation strengths. 

Figs. and capture all global features of the VT 
phase diagram. The detailed structure of the VT thresh¬ 
old manifold 7 pt(^O: 0) depends on the other two pa¬ 
rameters {X^P). We emphasize that the VT-symmetry 
breaking threshold is maximum when the distanee between 
the gain-site and the loss-site is maximum. Starting 
from gain-loss sites nearest to each other, mo ^ N/2, 
we expect that when the distance between them d = 
(A^ + 1 — 2mo) is increased, the T^T-symmetry breaking 
threshold 7 pt(<^) will decrease. Our results, however, 
predict otherwise. This surprising finding is due to open 
boundary conditions that ensure complete reflection at 
the two ends of the lattice. 

What is the origin of the positive threshold 7 pt when 
the underlying Hermitian Hamiltonian is not reflection 
symmetric? In the following section, we present a hidden 
symmetry of the eigenfunctions of H that is instrumental 
to a non-vanishing threshold. 
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FIG. 3. (Color online): 'PT-threshold 7PT(mo,0)/J dependence on the tunneling period p = 1//3 and modulation strength A. 
(a) 7 pt(^o, 0) for an N = 111 lattice with A = 0.5 shows a maximum when gain-loss sites are farthest from (mo = 4) or closest 
to (mo = 52) each other, (b) The same qualitative behavior is observed for the same lattice with a longer tunneling period 
1//3 = 7. (c) Results for an iV = 61 dimer lattice, /S = 1/2, show that as the modulation strength A increases, the threshold 
7pt is monotonically suppressed from its value in the A ^ 0 limit [48j . 


III. ORIGIN OF THE POSITIVE THRESHOLD: 

HIDDEN SYMMETRY OF THE AAH MODEL 

Let us recall how a positive T^T-breaking threshold 
arises in the case of traditional T^T-symmetric Hamil¬ 
tonians. If the tunneling Hamiltonian is T^T-symmetric, 
so are its eigenfunctions fa{k) with energies Cq,. In the 
presence of a VT symmetric potential, Eq.([^, the first- 
order perturbative correction to the eigenenergies is 
given by 

= iji\faimo)\^ - |/a(mo)7. (7) 

Since the eigenfunctions fa{k) have equal weights on 
reflection-symmetric sites (mo,mo), this correction as 
well as all higher odd-order corrections vanish for for 
all gain-locations, A^’^+^( 7 ,mo) = 0 for every mo [49] . 
This property ensures a real spectrum 60 ,( 7 ) poten¬ 
tial strength 7 < 7 pt- But what if Hq is not reflection- 
symmetric? Its arbitrary eigenstate fa{k) with energy 
satisfies the following difference equations at reflection- 
symmetric sites (k^k). 

tk-lfa{k - 1) + tkfa{k + 1) = -eafa{k), (8) 

tk-lfa{k - 1) + tlfa{k + 1) = -€afa{k), (9) 

where open boundary conditions are implemented by us¬ 
ing to = 0 = tTv- It follows that if eigen¬ 

functions, in general, will not have equal weights on the 
reflection-symmetric sites, |/o(^)| 7 ^ |/a(^)|- 

Figure shows typical eigenfunctions of the Hamilto¬ 
nian Ho{X, [3^ (j)) when the lattice size satisfies A + 1 = 0 
mod p. Panel (a) in Fig. shows the ground-state (GS) 
wave functions fcik) for an A = 29 site lattice with 
tunneling period p = 3 (solid line with open circles) and 




FIG. 4. (Golor online): Hidden symmetry of eigenfunctions of 
Ho{X, (j)). (a) A = 29 lattice with A = 0.5 shows reflection- 
asymmetric ground state (/3 = 1/3,0 = 0: solid line; /? = 
1/5,0 = tt: dashed line). However, the amplitudes on sites 
k = 0 mod (1/0) show reflection symmetry about the lattice 
center (red solid circles, blue solid squares), (b) Lowest energy 
midgap states for an A = 15 lattice with A = 0.5, 0 = 0, and 
0 = {1/2,1/4} show that their wave functions vanish at sites 
k oc 1/0. These states are unaffected by the balanced gain- 
loss potential. 

p = 5 (dashed line with open squares). The results are 
for tunneling modulation strength A = 0.5 and phases 
0 = 0, TT respectively. As is expected, both ground state 
profiles are reflection-asymmetric about the center site 
Tie = 13. However, these wave funetions have the follow¬ 
ing hidden symmetry. Solid red circles show the p = 3 GS 


















5 


amplitudes at sites k = 0 mod p = {3, 6,...}, whereas 
the solid blue squares show the p = 5 GS amplitudes at 
sites k = 0 mod p = {5, 10,..lu both cases the wave 
fuuctiou weights satisfy |/a(^o)| = |/a(^o)| if aud ouly if 
rriQ is au iuteger multiple of the tuuueliug modulatiou pe¬ 
riod p. This result is true for all eigeustates of i^o(A, /3, 0) 
if aud ouly if the lattice size N satisfies N 1 = 0 
mod p. It eusures that the reflectiou-couuterpart site 
iudex ffiQ = N ^ 1 — rriQ is also au iuteger multiple of 
the tuuueliug modulatiou period. This hiddeu symmetry 
is iustrumeutal to a positive VT threshold that we ob¬ 
serve wheu the gaiu poteutial is located at sites mo = 0 
mod p. It also implies, via the perturbatiou theory ar- 
gumeuts [49] , that the eigeufuuctious of the total Hamil- 
touiau H = Hq + V coutiuue to have this symmetry for 
7 ^ 7 pt- 

Next, we cousider implicatious of this hiddeu symme¬ 
try to localized midgap states that are, iu some cases, 
topological iu uature. Pauel (b) iu Fig.j^shows the lowest 
euergy midgap state for au A/" = 15 lattice with /3 = 1/2 
(solid hue) aud ^ = 1/4 (dashed hue); the results are 
for A = 0.5 aud phase 0 = 0. These states are localized 
at oue eud of the lattice. The surprisiug feature, shared 
by all localized midgap states, is the preseuce of uodes 
precisely at sites mo = 0 mod p = {p, 2p,...}. Wheu 
the lattice size satisfies N 1 = 0 mod p, due to the 
tuuueliug amplitude asymmetry at the two euds of the 
lattice, it follows that a midgap state must be localized 
at oue eud or the other, but uot equally at both euds. 
Therefore, the hiddeu symmetry discussed iu the previ¬ 
ous paragraph implies that its wave fuuctiou must vauish 
at sites mo = 0 mod p. This result is true for all (p — 1) 
localized midgap states. This remarkable property of the 
localized states shows that a balauced gaiu-loss poteutial 
will have uo effect ou them. Iu particular, the euergies 
of these states remaiu real aud these localized, topolog¬ 
ical states [41] remaiu robust eveu wheu the gaiu-loss 
streugth exceeds the threshold, 7 > 7 pt(A, 0, 0; mo). Iu 
receut years, the preseuce or abseuce of topological iu- 
sulator states iu T^T-symmetric Dirac aud SSH models 
has beeu exteusively studied [421145] . Our results show 
that robust, topological states occur iu a wide class of 
reflectiou-asymmetric Hamiltouiaus with a positive VT 
breakiug threshold. 


IV. ANALYTICAL AND BEAM PROPAGATION 
METHOD RESULTS: N=5 DIMER LATTICE 


The smallest experimeutal realizatiou of a lattice with 
reflectiou asymmetry aud a positive T^T-breakiug thresh¬ 
old, say iu a coupled waveguide array, will require N = b 
waveguides with a dimer tuuueliug profile ti = J (1 — 
Acos 0 ]) aud t 2 = J{1 + Acos 0 ), aud gaiu-loss poteutials 
at reflectiou-symmetric locatious (mo, mo). This au- 
alytically solvable case provides further iusights iuto the 
results preseuted iu this paper. 

It is easy to check that wheu mo = 1, the character¬ 


istic equatiou for the 5x5 Hamiltouiau H = Hq + V 
has complex coefhcieuts aud therefore, the VT threshold 
at mo = 1 is zero. Wheu mo = 2, the correspoudiug 
equatiou is giveu by 

X + 2 t\ — 7 ^) + (tf + ^2 + ^ 1 ^ 2 )] = 0 - ( 10 ) 


It follows from Eq.(lO) that the eigeuvalues of ff(j) are 
either purely real or occur iu complex coujugate pairs, 
aud the real spectrum has a particle-hole symmetry iu the 
T^T-symmetric phase, 7 < jpt [50]. The threshold gaiu- 
loss streugth at which the eigeuvalues trausitiou from real 
to complex coujugate pairs is giveu by 


7pt(A, 0) — J y^2(t^ + 


1-Jl- 


61^2 


(tlTtlY 


1/2 


( 11 ) 

It follows from Eq.(|ll|) that the threshold ypp is iuseusi- 
tive to the tuuueliug modulatiou A wheu ti ~ ^ 2 , aud it 
is maximally suppressed wheu cos 0 = ±1 (see Eig. ic)). 

It is also straightforward to show that the uuuormal- 
ized, zero-euergy, ’’edge-state” eigeuvector is giveu by 
I/) = (t|, 0 , — 0 , 1 )^ where = ^ 2/^1 is the ratio of 
tuuueliug amplitude ou the secoud boud to the tuuueliug 
amplitude ou the first boud m- Thus, the edge state has 
uodes at sites k = 0 mod 2 = {2,4}, is localized at the 
left eud (right eud) of the lattice wheu > 1 (t^ < 1 ), 
aud remaius uuaffected by the VT poteutial. 

Next, we test the validity of our predictious via beam- 
propagatiou method (BPM) [52] |53]. This method al¬ 
leviates the coustraiut of tight-biudiug approximatiou, 
by takiug iuto accouut the uouzero spatial dimeusiou 
of a ’’siugle site”. With a realistic waveguide-array iu 
miud [54], we obtaiu the time-evolutiou of au iuitially 
uormalized wave packet localized iu the ceuter waveg¬ 
uide, ^(x, 0) = exp [—{x — xs)‘^/4a‘^] /(27rcr^)^/^. Here x 
is the coutiuuous coordiuate trausverse to the waveguide 
array, x^ is the ceuter of the third waveguide, aud the 
iuitial wave packet size a is set to half-the-width of the 
waveguide, a = Wg/2. We remiud the Reader that iu 
this realizatiou, the wave fuuctiou 'ip{x,t) represeuts the 
slow-varyiug euvelope of the electric field E{x^zT) that 
also has a rapidly varyiug part proportioual to exp(//co 2 ;). 
The time-evolutiou of %jj{xT) is giveu by the Maxwell 
wave-equatiou iu the paraxial approximatiou mm. 


.90 




9^0 


dt 2 konQ dx‘^ 


+ cko 


n{xY 


^0 J 


0 . ( 12 ) 


Here c is the speed of light iu vacuum, no is the claddiug 
iudex of refractiou, n{x) = no + An{x) is the positiou- 
depeudeut iudex of refractiou, aud the iudex coutrast 
An ^ 10“^ 7 ^ 0 ouly withiu each waveguide. Eor 
An ^ no, the effective poteutial is proportioual to the 
iudex coutrast, V{x) oc An. We implemeut the gaiu aud 
loss by addiug imagiuary parts ±27 to the real iudex cou¬ 
trast An. Wheu the poteutial V(x) is uot real, the time 
evolutiou of the wave packet is uot uuitary, aud there¬ 
fore the total iuteusity I{z) = /dx| 0 (x,t = zno/c)p is 
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FIG. 5. (Color Online) BPM results for an array of iV = 5 coupled waveguides. The width of each waveguide is Wg = 5 //m and 
the distances d ~ 10 /xm between the waveguides are chosen such that the tunneling ratio satishes ti/t 2 = 3/5. The left-hand 
panel in each row shows the space-and-time dependent intensity I{x,z); the initial state is a normalized Gaussian centered on 
the 3rd waveguide (shown by a white half-circle). The right-hand panel in each row shows the net intensity I{z) as a function 
of time t, or equivalently, the distance = ctino along the waveguide. The center panel shows the schematic of gain and loss 
locations, (a) The array is in the VT broken phase when the gain-location is first waveguide, (b) With the same gain-loss 
strength, the array is in the PT-symmetric phase when mo = 2. (c) When the gain is doubled, the system transitions into the 
VT broken phase. 


not a constant; note that we have switched to the dis¬ 
tance along the waveguide 2 ; = ct/no as a stand-in for 
the time, for it allows an easier comparison with typical 
experimental setups. 

Figure shows the results of such an analysis. Each 
row shows the space- and time-dependent intensity 
I{x,z) = |'0(x,^)p (left-hand panel); a schematic of the 
corresponding 5-site lattice (center panel); and the time- 


dependence of the total intensity I{z) (right-hand panel). 
The first row shows that when the gain and loss potentials 
Tiy are located on the first and the last waveguides re¬ 
spectively, the net intensity I{z) increases monotonically 
with time, indicating a T^T-symmetry broken phase. The 
second row shows the results when the gain-loss are at 
sites mo = 2 and mo = 4 respectively. It is clear from the 
I{x^z) plot that the wave packet undergoes oscillations 
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across the lattice along with some amplification. This pe¬ 
riodic behavior is also manifest in the total intensity I{z) 
and shows that the system is in the T^T-symmetric phase 
for the same gain-loss strength. The bottom row shows 
that when the gain-loss strength is doubled, the system 
enters T^T-broken phase, as evidenced by monotonically 
increasing net intensity I{z). 

These results demonstrate that the non-trivial depen¬ 
dence of the T^T-breaking threshold on the gain location 
mo for reflection-asymmetric models is robust. 


V. DISCUSSION 

In this paper, we have discovered that a broad class 
of Aubry-Andre Harper models isain] with reflection- 
asymmetric tunneling profile can have a positive VT- 
symmetry breaking threshold. This occurs when P = q/p 
is rational, and the lattice size N and the gain-potential 
location mo < N/2 both satisfy + 1 = 0 mod p 
and mo = 0 mod p. These constraints ensure that the 
loss-potential location m also satisfies m = 0 mod p. 
Through the tight-binding analysis of the lattice model 
and a BPM analysis of its continuum counterpart, we 
have shown that our predictions remain valid for realistic 


waveguide arrays. The A AH lattice models investigated 
here are known to support topological states IH]. They 
thus provide an avenue to experiments in which the in¬ 
terplay between T^T-symmetry breaking and topological 
properties can be studied. 

We note that this paper is based on an effective, single¬ 
particle Hamiltonian that permits amplification and de¬ 
cay. Prima facie^ these results predict the existence of 
topological insulators with positive VT breaking thresh¬ 
old [42l [43] , since our model makes no reference to the 
quantum statistics of the particle. In reality, however, 
amplification of a single degree of freedom is incompati¬ 
ble with the Pauli principle. Thus, our results can apply 
to fermions only if the gain and loss are associated with 
the bulk Fermi sea, and not with a single quantum de¬ 
gree of freedom. In the bosonic case, amplification of a 
single degree of freedom is permitted and our results are 
directly applicable. We have also ignored two-body inter¬ 
actions; they become important only in the T^T-broken 
phase as the on-site intensity (light) or density (massive 
bosons) is amplified. 
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